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Polynomial convexity of real submanifolds of C" is a well-studied subject in complex analysis 
due to its deep relation to the approximation problems, pluripotential theory and Banach algebras 
(see, for instance, |2, 26] for a detailed discussion). M. Gromov |15) found remarkable connections 
between the polynomial (or the holomorphic disc) convexity of real manifolds and global rigidity 
of symplectic structures. In the present work we prove that a generic Lagrangian surface in 
is polynomially convex near an isolated singularity which is topologically an unfolded Whitney 
umbrella. This study is inspired by the work of A. Givental |14] where he proved that a wide 
class of compact Lagrangian surfaces in a symplectic 4-manifold can be realized as a compact 
surface in the complex affine space equipped with the standard symplectic form. This surface 
is Lagrangian outside a finite number of complex points (i.e., points singular for the CR structure, 
see below) and unfolded Whitney umbrellas. 

Denote hy z = x + iy and w = u -^riv the standard coordinates in . Let 



be the standard symplectic form on C^. A smooth map : — ?• is called symplectic if 
(f)*uj = oj. If a symplectic map is a (local) diffeomorphism, we call it a (local) symplectomorphism. 
A smooth map i : 5" — )• (C^,ti;) from a smooth real surface S is called isotropic if l*uj = 0. 
A. Givental |14) showed that near a generic point p G S, which is an isolated singular point of l 
of rank one, the map 



1. Introduction 



uj = dx f\ dy + du f\ dv 




(1) 
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is a local normal form for i. In particular, this means that there exists a local symplectomorphism 
near sending l{S) onto a neighbourhood of the origin in S := 7r(M^). The set S, as well as 
i{S) near is called the unfolded (or open) Whitney umbrella. Our main result is the following. 

Theorem 1. Suppose : — )• is either a generic real analytic symplectomorphism near the 
origin, or the identity map. Then there exists a neighbourhood of the point </)(0) in the surface 
with the compact polynomially convex closure. 

The case where (p is the identity map is considered separately since it is not generic. This 
implies that the Whitney umbrella S is polynomially convex near the origin. The above theorem 
also holds under weaker assumptions, namely, if </> is a generic local real analytic diffeomorphism 
and D(j){0), the differential of (j) at zero, is symplectic, or if (/> is a C°°-smooth symplectomorphism 
with the jet at the origin satisfying some additional assumptions. See Section 5 for details. 

Denote by M[p,r) the open Euclidean ball of C"^ centred at p and of radius r > 0. As an 
application of Theorem [1] we obtain the following result. 

Corollary 1. Let (j) be as in Theorem d Then for e > sufficiently small, any continuous 
function on nB(</>(0),e) can be uniformly approximated by holomorphic polynomials. 

For a real n-dimensional submanifold E (closed or with boundary) in C" the polynomial con- 
vexity of E is connected with its global topological properties. Generically E is totally real, i.e., 
at a generic point the complex span of its tangent space coincides with C" . It is well-known [21 [26] 
that an n-dimensional totally real submanifold of C" is locally polynomially convex, i.e., for every 
point p £ E there exists r > such that the intersection E n M{p,r) is polynomially convex. 
This is no longer true in the global situation. Indeed, M. Gromov p5] proved that a compact 
Lagrangian (with respect to the standard symplectic structure on C") submanifold of C" contains 
a boundary of a non-constant holomorphic disc. Extending his ideas, H. Alexander [Ij proved 
that for every totally real submanifold E of C" of real dimension n there exists a non-constant 
bounded holomorphic disc with the boundary attached to E except possibly a single point of the 
unit circle. However, the methods used to prove these results cannot be immediately generalized 
to the case when E has singularities. An interaction between the rational convexity and the 
symplectic structure was studied by J. Duval ^lOj. Recently J. Duval and D. Gayet [12] studied 
the rational convexity of certain immersed Lagrangian manifolds. 

The local polynomial convexity can fail near a point where E is not totally real. Generically 
these points are isolated in E. In the complex dimension n = 2, the tangent space of is a 
complex line, so such points are called complex. There are three types of complex points: elliptic, 
hyperbolic and parabolic (see, for instance, [21 [26]), and the local polynomial convexity depends on 
the type. H. Bishop f5] and C. Kenig - S. Webster p9] proved that a neighbourhood of an elliptic 
point in E has a non-trivial hull. On the other hand, F. Forstneric and E. L. Stout [13j proved 
that E is locally polynomially convex near a hyperbolic point. The parabolic case is intermediate 
and in general both possibilities occur. This case was studied by B. Joricke [161117]. Normal forms 
of real analytic totally real or Lagrangian surfaces near complex points are studied by J. Moser - 
S. Webster pQj and S. Webster p7]. The present work shows that the local convexity properties 
near a Whitney umbrella are similar to the case of a hyperbolic point. 

We thank S. Nemirovskii and V. Shevchishin for bringing our attention to this problem and for 
helpful conversations. 
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2. Geometry of Whitney umbrellas 

The map vr : R?, , — >• Mf , ^ given by (fTl) is a smooth homeomorphism onto its image, non- 
degenerate except at the origin, where the rank of vr equals one. It satisfies it*uj = 0, and so S is 
a Lagrangian submanifold of (C'^^uj) with an isolated singular point at the origin. Thus, 

S = {{z, w) ^ : X = ts, u= -—, y = t"^, v = s; t, s G M}. 

The crucial role in our approach is played by an auxiliary real hypersurface M defined by 

M = {(z, u-) G : p{z, u>) = - yv^ + - y^ = 0}. (2) 

Clearly, S is contained in M. Note that the hypersurface M is smooth away from the origin, and 
strictly pseudoconvex in B(0,e) \ {0} for e sufficiently small. 

Suppose now that (p-.C^^C^ is a local smooth diffeomorphism near the origin such that its 
linear part D(f){0) at the origin is a symplectic map. Without loss of generality we may assume 
that (/)(0) =0. The standard symplectic structure on is given by the matrix 

n=( 

where I2 denotes the identity matrix on M?. Similarly, we write 

Dm = l^c d)- 

The condition that D(f){0) is symplectic means that (D(j){0)Y 0, D(j){0) = 0, (where t stands for 
matrix transposition). Therefore, the real (2 x 2)-matrices A, B = (bjk), C, D = [djk) satisfy 

A^D - &B = I2, A*C = &A, D^B = B*D. (4) 

The standard complex structure of in real coordinates is given by the matrix 



J 



-h 
h 



which corresponds to multiplication by i. We perform an additional complex linear change of 
coordinates ip. Let -i/; : — )• be a linear transformation given by the 4x4 matrix 

-B* 
5* D 



t nt I ■ (5) 



This matrix commutes with J and so gives rise to a non-degenerate complex linear map in C^. 
Let T,' = ip o 0(5]) and M' = ip o (f){M). The differential at the origin of the composition -0 o is 
given by 

D{^ocPm=(^^^ (6) 

where we used identities ^ to simplify the matrix. Further, a direct calculation shows that 

G = {gkj) = B'B + D'D, (7) 

and therefore, the matrix G is symmetric with positive entries in the main diagonal. The deter- 
minant 

A = 511522 - 5i2 (8) 
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of G coincides with that of the matrix in ([5]) corresponding to a C-Hnear map of €?. Hence A is 
also positive. Let p' = p o {ip o , and 

= {(z', w') G : p'{z', w') < 0}. (9) 

It follows from © and ([6]) that 

p'iz',w') = x'^ + ^n'' + o(|(z',u;')P). (10) 

In particular, the function p' is strictly plurisubharmonic in a neighbourhood of the origin, and 
the hypersurface M' is strictly pseudoconvex in a punctured neighbourhood of the origin. 

Lemma 1. The polynomial hull of the set 18(0, e) n S' for sufficiently small e > is contained in 
O'nB(0,e). 

Proof. Choose e > small enough such that p' is strictly plurisubharmonic in B(0,e). The 
polynomially convex hull of 18(0, e) n S' is contained in 18(0, e). By a classical result (see, for 
instance, |18]). the polynomially convex hull of 18(0, e) n S' coincides with its hull with respect to 
the family of functions plurisubharmonic in B(0,e). Since for any point p in B(0,e) \ 17 , we have 
p'{p) > 0, the assertion of the lemma follows. □ 

3. Characteristic foliation and polynomial convexity 

Let X be a totally real surface embedded into a real hypersurface y in C^. Define on X a field 
of lines determined at every p £ X hy 

Lp — TpX n HpY^^ 

where HpY = TpY r\i{TpY) denotes the complex tangent line to Y at the point p. Integral curves, 
i.e., curves which are tangent to Lp at each point p, of this line field define a foliation on X. It is 
called the characteristic foliation of X. 

We consider the characteristic foliation of S\{0} C M and {tjjo(f)){Ti)\{{f\ c 4>{M). Characteris- 
tic foliations are invariant under biholomorphisms. Therefore, in order to study the characteristic 
foliation on with respect to (j){M), it is sufficient to study the characteristic foliation of 

Y,' = lp o <j)(Ti) with respect to M' . Considering its pull-back by o o ^i" we obtain a smooth 
vector field in a neighbourhood of the origin in M^^ with the singularity at the origin. For 
simplicity, its trajectories will also be called the leaves of the characteristic foliation. 

Our ultimate goal is to prove the following 

Proposition 1. Let (j) he as in TheoremUl There exist two rectifiahle curves 71 and 72 passing 
through the origin in M^^ with the following property: if K C M^^ is a sufficiently small compact 
containing the origin and not contained in 71 U72, then there exists a leaf j of the characteristic 
foliation on S' such that K j ^ but K does not meet both sides of 'j. 

The proof of the proposition will be given in Sections H] - [71 It is based on the local theory of 
dynamical systems and can be read independently from the rest of the paper. 

Assuming Proposition [1] we now prove our main results. The proof is based on the argument due 
to J. Duval [llj and B. Joricke [16l[l7j. Suppose that (p satisfies the assumptions of Theorem [H 
and S' = (/"(S). First we establish non-existence of holomorphic discs attached to T,' near the 
Whitney umbrella. By a holomorphic disc we mean a map / : A — t- holomorphic in the unit 
disc A C C and continuous on A. As usual, by its boundary we mean the restriction /|aAj we 
identify it with its image. 
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Corollary 2. There exists e > with the following property: a holomorphic disc / : A — )• B(0,e) 
with the boundary attached to S', i.e., satisfying f{dA) C S', is constant. 

Proof. Fix e > such that the function p' in (|10p is strictly plurisubharmonic in the ball ]B(0, 2e). 
Suppose that / is not constant. Then the function p' o f \s subharmonic in the unit disc and 
the maximum principle implies that /(A) is contained in {p' < 0}. By the uniqueness theorem 
the set of points /^^(O) has measure zero on the unit circle. Since E' is totally real outside the 
origin, it follows by the boundary regularity theorem [7] that / is smooth (even real analytic) up 
to the boundary outside the pull-back /~^(0). Applying the Hopf lemma (see for instance [23 1) to 
the subharmonic function /o' o / on A we conclude that / is transverse to the hypersurface M' at 
every point different from the origin. Therefore, the complex line tangent to /(A) at a boundary 
point is transverse to the tangent complex line of M' at this point. In particular, the boundary 
K := f{dA) is transverse to the leaves of the characteristic foliation of S'. This contradicts 
Proposition [1] □ 

Proof of Theorem [IJ Given a compact set K, we denote by K its polynomially convex hull. Let 

X = S'n]B(0,e), x'''" = x\x, x*^=X^nx. 

A local maximum principle of Rossi [24J states that if ii' is a compact set in C", E C K is compact, 
U is an open subset of C" that contains E, and if / G 0{U), then = ||/||(_Bn_K')u9£;) where 

the boundary of E is taken with respect to K. By choosing E = X^'^^ and {/ = we see that 
X^'^ = X^'^*. Therefore, to prove that X is polynomially convex, it is enough to show that is 
polynomially convex. 

By an analytic curve in an open set U in we mean an irreducible complex 1-dimensional 
analytic subset in U . 

Lemma 2. [Duval, [llj/ Let p £ X \ {0}, and let Q' be given as in ([9]). Suppose there exist two 
continuous families {Vt}t&[o,i) o-^d {Wt}te[o,i) of analytic curves in an open neighbourhood of Q' 
with the following properties: 

(i) Vq and Wq meet X transversely at p and with opposite signs of intersection, and 

(ii) for t > 0, the varieties Vt and Wt are disjoint from X*''. 
Then pi X^\ 

Duval's original result is stated for the 0(r2)-hull of a smooth totally real surface X C 90, 
where 57 C is a strictly pseudoconvex domain. However, the proof is also valid in our situation. 
Indeed, suppose p £ X \ {0}, and {Vj}, {Wt} are as in Lemma [2j Duval's construction produces 
an arbitrarily small neighbourhood U of p and a complex 1-parameter family {Cq} of analytic 
curves in U such that Ca fill out U\X, avoiding X, and such that each Ca can be swept out of Q' 
through a continuous family. Using this family we are now in the position to prove that p ^ X*^'. 
For that we apply the classical characterization of polynomially convex hulls due to Oka (see |21j . 
|25j . |26j). which we state in the form suitable for our purposes: 

For a compact set X C C'^, a point p E does not lie in the polynomially convex hull X of 
X, if there exists a neighbourhood U of p and a continuous family {St}te[o,i] of analytic curves in 
U with the properties that 

(i) p G ^O; 

(ii) St are disjoint from some open neighbourhood of X, 

(iii) for any compact K C U , there exists tx < ^ such that St O K = for tx < t < 1, and 

(iv) there exists to < 1 such that {St} (1 X = for to < t < 1. 
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In our situation we take any point in U \ X and choose the family {St} to be a continuous 
subfamily of {Ca} that initially passes through this point and then leaves fi'. Note that by 
Lemma [H the polynomial hull of X is contained in 0', and therefore condition (iv) in Oka's 
characterization holds. This shows that no point near p can be in X'^^^, and therefore p does not 
belong to X*''. This verifies Lemma [21 

Let 7i and 72 be as in Proposition [H and suppose that X*''" is not contained in the union of these 
curves. Then, by Proposition [1], there exists a leaf 7 of the characteristic foliation that touches 
X^^ at a point p G X \ {0}. There exists a small neighbourhood U of p and local holomorphic 
coordinates centred at p such that in these coordinates M n ?7 is strictly convex. Again, the local 
analysis of Duval [Tl] produces two families of analytic curves in U that satisfy the conditions of 
Lemma El The intersection of each curve with is compactly contained in U, and so they can 
be considered as curves in a neighbourhood of O'. By Lemma O p ^ X^^'. 

It follows from the above considerations that X*'' is contained in the union 71 U 72 . A rectifiable 
arc is polynomially convex; moreover, if Y is compact and polynomially convex, and F is a compact 
curve, then the set (y U P) \ (y UP) is either empty or contains a purely one dimensional analytic 
subvariety of \ (y U P) (see [26], p. 122). By taking Y and P to be our rectifiable curves jj, we 
see that their union cannot bound a one dimensional variety, and therefore X*^ is polynomially 
convex. Hence, X^^'^ is empty, which implies that X is polynomially convex. Theorem [1] is 
proved. □ 

Proof of CorollarylJl Let (f){0) = p. By Theorem [1] there exists e > such that 0(S) n M{p,e) is 
polynomially convex. We may further assume that n dM{p, e) is a rectifiable curve. By the 
result of J. Anderson, A. Izzo, and J. Wermer [3^ Thm. 1.5], if X is a polynomially convex compact 
subset of C", and Xq is a compact subset of X such that X \ Xq is a totally real submanifold of 
C", of class C^, then continuous functions on X can be approximated by polynomials if and only 
if this can be done on Xq. We apply this result to X = (/>(S) n M{p,e) and Xq = {p} U ((/)(S) n 
dM(p,£)). Indeed, Xq, is polynomially convex, and furthermore, by [26j, Thm. 3.1.1 and Cor. 
3.1.2, continuous functions on Xq can be approximated by polynomials. From this the corollary 
follows. □ 

The rest of the paper is devoted to the proof of Proposition [H 



In this section we deduce the dynamical systems describing the characteristic foliations on S 
and E'. In Sections [6] and [7] we will discuss the topological behaviour of these foliations near the 
origin. 

4.1. Foliation on E. The tangent plane to S \ {0} is spanned by the vectors 



4. Reduction to a dynamical system 




The directional vector of the characteristic line field is determined from the equation 



X = aXt + /3X 



(11) 
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where a = a{t, s),f3 = f3{t, s) are some smooth functions on M? \ {0}, and the vector X belongs 
to the complex tangent H^(j. g^M. Let 

(J ?)' '^={1 ~Q 

Multiplication by i of a vector in corresponds to multiplication by J of the corresponding 
vector in M^. For v G TpM, the inclusion v G HpM holds if and only if v,iv G TpM. Therefore, 

X e ^ {J{aXt + pXs),Vp) = 0, 

where (•,•) is the standard Euclidean product in M^, and Vp is the gradient of the function p. 
Therefore, we can choose 

a={JXs,Vp), P = -{JXt,Vp). (12) 

A calculation yields 
and 



Vp = {2ts, 3^^ -s^ - 3i^, -2t'^s), 

a = -3t^ -ts^ - 3^^ 
P = s^ + 4fs + 7st^. 



Thus, 



(14) 



X = aXt + pXs = adTT (^^^^ + PdTT (^^^^ = dTT (^"^^ , (13) 

where dvr is the differential of the map vr. It follows that the characteristic foliation on S \ {0} 
(or, more precisely, its pull-back on M^\{0} by the parametrization map vr) is given by the system 
of ODEs of the form 

i = -3*3 - ts'^ - 3t^ 
s = + 4*^5 + 7st 
where the dot denotes the derivative with respect to the time variable r. 

4.2. Foliation on S'. Let / : M"^ be given by 

f := tp o (j) o TT, 

where we use the notation of the previous section. The directional vector of the characteristic 
foliation on S' is determined by 

X' = aX't + pX'^, 

where = df/dt and X'g = df/ds, and a = a{t,s),f3 = f3{t,s) are some smooth functions 
on \ {0} which are chosen in such a way that the vector X' belongs to the complex tangent 
Hf{f^s-)M'. We have 

X' e Hf^t,s)M' ^ {J{aX[ + Vp') = 0, 
where p' is a defining function of M' , and the gradient Vp' is expressed in terms of (i, s) using 
the parametrization /. Therefore, we can choose 

a{t,s) = {JX'„Vp'), P{t,s) = ~{JXt,Vp'). (15) 

Thus, 

X' = aXl + f3Xi = df(^''^y (16) 

It follows that the characteristic foliation on S' is determined by the system of ODEs of the form 

i = a{t,s) . . 
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We write f{t, s) = {fi{t, s), . . . , /^{t, s)), where using Q and ([T]) we may express each fj as a 
power series in (t, s): 

f^{t,s) = x+ lU^xH^y'v"^ = 

j+k+l+m>2 

ts + fl2S^ + fUs^ + fit's + fl^s^+ fjkt's", (18) 



j+k>4: 



where fjj^i,^ and fjj^ are real numbers. Similarly, 



j+k+l+m>2 

^t^ + /oV + /iW + /2V^^+/oV+ E fh^'^'- (19) 

j+k>4 

Denote by Cj^ the entries of the matrix E in Then 

hit, s) = eux + ei2U + guy + guv + ^ ff^i^ xH^y^v^ = 



i+k+l+m>2 

2ei_2. 
3 



j+fc>4 



/4(t, s) = e2ix + e22U + guy + g22V + 



EnA j k I m 

fjklmX^U yv 



j+k+l+m>2 
2622 _ 

3 



522S + 5i2i' + e2it5 + /oV + ^i' + /iW + /2V'^ + /oV+ E ffk*'''- (21) 

j+k>4 



From these formulas we immediately obtain 

/ s + 2flts + fl^s' 

* 2giit + ens + 2612*^ + 2fl^ts + ff^s'^ 
\ 2gut + 6213 + 2622*2 + 2/lits + ff^s'^ ) 



and 



/ t + 2fl^s + fl^e + 2/1^5 + 3/I3S2 
2/oV + /2V' + 2/iV^ + 3/oV 

512 + exit + 2/32S + + 2fl^tS + S/^gS^ 

V 522 + e2lt + 2/^28 + + 2/42ts + 3/4352 / 

-1 



(22) 



(23) 



The defining equation of M' can be chosen to be /f? o (^ o (^) ^, where /O defines M as in ([2]) 
Let (x', n', y', v') be the coordinates in the target domain ip o (p, in particular, we have x' = fi 



u' = f2, y' = /a, and v' = f^. Let 

{Di^P o ^)(0))-i = ( ^2 0, ^ , i^' = (e;.,), G' = (5;,). 



(24) 
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Then 

y j+k+l+m>2 j+k+l+m>2 

e'nx' + e'i2n' + ^iiy' + 9'i2v' + J2 ^%im x'H'^y'^v'"', 

j+k+l+m>2 

e',,x' + e',,u' + g[,y' + g',,v' + ^%im ) . (25) 

j+k+l+m>2 J 

Therefore, 

p'ix', u', y', v')={x'+ Y ^Um x'^u'Vv'^ 

j+k+l+m>2 



2 



2 

3 



+ Zh'+ E ^\z„^^VVV"^| -(elix' + .-. + ^^' + E^^^n^^VVV'")" (26) 



4 

\ j+k+l+m>2 I 

Note that in (j26p the only quadratic terms are x'^ and By taking partial derivatives in the 

above expression with respect to x' ^u' ,y' and v' , and expressing the resulting vector in terms of 
{t, s) we will obtain the coordinates of the vector 

To determine the phase portrait of the characteristic foliation we will only need some low order 
terms in the power series 

a{t,s) = E (^jkt^s^, (i{t,s) = E l^jkt^^'. 

j,k>0 j,k>0 



Therefore, instead of explicit differentiation of (j26p . we will employ a different strategy for com- 
puting coefficients of the terms of lower degree in the {t, s)-Taylor expansion of a and (3. 

4.3. The power series of a. We have 

a{t,s) = {JX'^Vp') = -{X'Js • Rx - {X'^)i • Ru + {Xi)i ■ Ry + {X'j2 ■ Rv (27) 

We proceed in several steps computing the coefficients in the expansion for a. To begin with, 
there cannot be a free term in the power series of a because every term in Vp' will necessarily 
have positive degree in i or s. 

Term t: Since no component of Vp' can contain a degree zero term or the monomial i, there is 
no term t in a. 

Term s: The first two components of X'^ do not contain free terms, therefore, monomial s can 
appear in a only if Rx or R^ will contain it. By inspection of (jlSp - (|2ip we see that y' and v' are 
the only terms that can produce monomial s. Therefore, for s to appear in R^ or R^, the function 
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p' must contain at least one of the terms x'y', x'v\ u'y' or u'v' . However, from (j26p neither of 
these terms exists. Thus, there is no monomial s in the power series of a. 

Term ts: We inspect terms in X'^ of degree lower than ts. These appear in {Xg)i (terms t and s), 
in {Xg)2 (term s), in {X'g)^ (a free term, t and s), and in (X^)4 (a free term, t and s). Therefore, 
for ts to appear in a, at least one of the following options must occur: 

(1) either or Ru has t, s or ts; 

(2) Ry has either t or s 

(3) Ry has t. 

Of the above three options only (1) can happen: p' contains the term x'^, and therefore, Rx 
contains 2ts. It follows now from (jl8p . (|23p and (j27p that an = —2gi2. 

To simplify further considerations, we note that term t cannot occur in any of the components 
of the vector V/o'. 

Term t^: By inspection of X'^, we conclude that either R^ or i?„ has term t^, so V/>' must have 
either x'y', x'v' , u'y' or u'v' , neither of which appears. This means that a does not contain term t"^. 
Term s^: By inspection of X'^, the following options are possible: 

(1) either Rx or R^ has s or s^; 

(2) either Ry or R^ has term s. 

Option (2) is impossible, but p' can have terms u''^, u'v'"^ or u'y''^ which gives (1). We have 
the following expression for ao2, which depends on the coefficients of the Taylor expansion for 

(V^o0)-l: 

"02 = 2(^0002922 + /02 + ^00205l2)- 

Term t^: By inspection of X'^, the following options are possible: 

(1) either Rx or Ru has at least one of or t^; 

(2) Ry has t^. 

Option (2) can happen only if p' would have y'^ or y'v', which is impossible. For the same reason 
in option (1) terms Rx or Ru cannot produce t"^. The only term in Vp' that can produce t^ is u'. 
Therefore, the only possibility in (1) is the term in Ru, which indeed happens since p' contains 
u'^. It follows that aso = —3522- 
Thus, 

a(t, s) = -2gi2ts + ao2S^ - 3^22*^ + ^ ajkth^ . 

j+k>2, {i,fc)^(3,0) 

4.4. The power series of /3. We have 

/3(t, s) = -{JX't, Vp') = (X;)3 • Rx + {X't)4 ■ Ru - {X't)i ■ Ry - (X^a • R^. 

Again, there cannot be a free term in /? because every term in Vp' will necessarily have positive 
degree in t or s. Further, no component in V/o' can produce a term t, and so the power series of 
/3 cannot contain monomial t. 

Term s: Since no component of X^ contains a free term, /3 cannot have monomial s. 

Term ts: By inspection of Xj- we conclude that either Rx or Ru must have term s, which is 
impossible. Hence, /3 does not contain monomial ts. 

Terms t^ and s^: Analogous considerations show that these terms cannot appear in /3. 
Term t'^s: By inspection of X^ the following is possible for R: 
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(1) Rx has at least one of t^, s, or ts; 

(2) Ru has at least one of t^, s, or ts; 

(3) Ry has t^] 

(4) Ry has s. 

Options (3) and (4) imply that p' has f'^, y'^, or u'y', neither of which is possible. Option (2) 
implies that p' has u'y' , u'v' and u'x'. Neither of these terms are present in p' , so (2) is also not 
possible. Option (1) implies that p' has at least one of x'y', x'v', or x''^. Only the latter happens, 
and so /32i = 4(jrii. 

Term ts^: This term can appear in /?. We have 

/3i2 = 2eii + 6512/02- 

Term t^: By inspection of X^, the only option is that either Rx or i?^ has term i^. This is however 
not possible. 

Term t^: The possibilities for R are as follows: 

(1) Rx has at least one of or t^; 

(2) has at least one of i^, or t^; 

(3) ii^ has t^. 

Option (3) cannot occur. The only possible option in (1) or (2) is that appears in Ru- This 
comes from the term u'^ in p' . It follows that /3o4 = ^912- 

Term s^: We have 

1 9 2 
Poz = 2eii/o2 + 2^2i/o2- 

Combining everything together we get 

s) = 4511^25 + /3i2ts2 + /5o3s3 + 6ffi2t^ + Yl Pjkth^. 

j+k>3, {j,k)=i{4,0) 



We note that if (p is merely a smooth diffeomorphism, then the above calculations give the values 
for the jets of a and (3 at the origin of the corresponding orders. In either case the characteristic 
foliation on S' is given by 

r i = a{t, s) = -2gi2ts + 0025^ - 3^22*^ + o{\t\^ + |sp + \ts\) , . 

\ s = (3{t, s) = Agnth + p^ts^ + (3o3S^ + 6512*^ + o{\t^s\ + \ts^\ + \s\^ + \t\^). 

It is easy to see that for a generic symplectomorphism (j) : {x,u,y,v) 1— )• {x' ,u' ,y' ,v') and a 
generic V the coefficients 002, /3i2) Poz do not vanish. Indeed, if il) is close to the identity map 
and the component u' of (f) contains the term av"^ with a 7^ 0, then 7^ and ao2, 1^12, P03 do 
not vanish. Therefore, they do not vanish generically. 

Remark. It follows from the above considerations that our restriction on (j) to be generic 
involves only the 2-jet of (f) at the origin. In other words, it suffices to require in Theorem [1] that 
(j) has a generic 2-jet at the origin. 

Lemma 3. Let cj) be a local symplectomorphism near the origin, and let X he the vector field 
near the origin in corresponding to the characteristic foliation on E'. Then X does not vanish 
outside the origin. 
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Proof. Since (j) is symplectic, 0(S\{O}) is a Lagrangian surface, in particular, totally real. There- 
fore, ip o (j)(T, \ {0}) does not contain complex points. Further, it easily follows from (jlSp that 
ct(io)*o) = /^(^OjSo) = 0, {to, So) 7^ if and only if f{to,so) is a complex point of S'. From this 
the result follows. □ 

5. Generalities on planar vector fields 

For the proof of Proposition [T] we need to determine the topological structure of the orbits or 
maximal integral curves associated with the vector fields defined by ()14p and ()28p . Both systems 
have higher order degeneracy (the linear part vanishes) at the origin, and consequently it is a 
nonelementary singularity of ()14p and l\28h . Therefore, standard results, such as the Hartman- 
Grobman theorem, do not apply here. Instead, we will use some more sophisticated tools from 
dynamical systems. We will be primarily interested in understanding the topological picture 
of ()14p and (128p near the origin up to a homeomorphism preserving the orbits. In this section we 
outline relevant results and recall some common terminology. 

The local phase portrait of a vector field near a nonelementary isolated singularity can be 
determined through a finite sectorial decomposition. This means that a neighbourhood of the 
singularity is divided into a finite number of sectors with certain orbit behaviour in each sector. If 
the vector field has at least one characteristic orbit (i.e., orbits approaching in positive or negative 
time the singularity with a well-defined slope limit), then the boundaries of the sectors can be 
chosen to be characteristic orbits. The overall portrait is then understood by gluing together the 
topological picture in each sector. The general result due to Dumortier [8] (see also [9j) can be 
stated as follows: 

Suppose that a -smooth vector field X singular at the origin in satisfies Lojasiewicz's 
inequality 

\X{x)\ > c|j;|^, c> 0, /c G N, 

for X G is some neighbourhood of the origin. Then X has the finite sectorial decomposition 
property, that is, the origin is either a centre (all orbits are periodic), a focus/node (all orbits 
terminate at the origin in positive or negative time), or there exists a finite number of charac- 
teristic orbits which bound sectors with a well-defined orbit behaviour (hyperbolic, parabolic, or 
elliptic). If the vector field X has a characteristic orbit, then its phase portrait is determined by 
its jet of finite order k, in the sense that any other vector field with the same jet of order k at the 
origin has the phase portrait homeomorphic to that of X . Further, whether the vector field X has 
a characteristic orbit depends only on a jet of X of some finite order. 

The original proof of the above result in ^ is based on the desingularization by means of 
successive (homogeneous) blow-ups. After each blow-up the singularity is replaced by a circle, 
and after a finite number of such blow-ups one obtains a vector field with only non-degenerate 
singularities. The construction of the blow-up maps depends only on a finite order jet of the 
original vector field at the origin. From the configuration of the singularities of the modified 
system on the preimage of the origin under the composition of blow-ups, it is always possible 
to deduce if the original vector field has a characteristic orbit. If such an orbit exists, then the 
singularity is not a centre or a focus, and the phase portrait is determined by a jet of finite order. 
Further, the Lojasiewicz inequality holds for any real analytic vector field in a neighbourhood of 
an isolated singularity (see, e.g., |3]) and, in particular, in our case, in view of Lemma [3l 

Alternatively, it is possible to use quasihomogeneous blow-ups, which are chosen according to 
the Newton diagram associated with X (see [22]). The advantage is that this gives a computa- 
tional algorithm for constructing the sectorial decomposition for a particular system. A detailed 
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discussion of this approach for real analytic systems is given in Bruno [6] in the language of nor- 
mal forms. Using Bruno's method we will show that for a real analytic 4> in general position, the 
vector field defined by (j28p will always have a characteristic orbit, and its phase portrait near the 
origin is a saddle. 

If in Theorem [1] the map (p is smooth, then the vector field corresponding to the characteristic 
foliation is only smooth, and the Lojasiewicz inequality imposes additional assumption on the 
vector field, and therefore on (p. The Lojasiewicz condition depends on the jet of the vector field 
at the origin and holds for all jets outside a set of infinite codimension in the space of jets, but it 
is not clear whether for a generic smooth symplectomorphism the inequality is satisfied. However, 
assuming that the Lojasiewicz condition does hold, the topological picture of the characteristic 
foliation is determined by its finite jet at the origin. Therefore, we may consider a polynomial 
vector field obtained by truncation of (j28p at sufficiently high order without distorting the phase 
portrait of the system. After that we may apply Bruno's method to determine its geometry. Thus, 
in Theorem [1] we may assume that is a generic smooth symplectomorphism such that the vector 
field corresponding to the characteristic foliation satisfies the Lojasiewicz inequality. 

If in Theorem [T] the map (/> is a real analytic diffeomorphism with D(j){0) symplectic, then all 
of the arguments go through provided that the vector field ()28p vanishes at the origin only. The 
latter holds for the following reason: consider near the origin the complexification F of the real 
analytic map f = tp o (j) o n : ^ C"^ . Then F : — )• is a holomorphic map such that 
F\^2 = f, in particular, F(R^) = S'. Moreover, since / has rank 2 outside the origin, it follows 

that the Jacobian of F does not vanish on M? \ {0}, and therefore, F is a local biholomorphism 
near any point on M? \ {0}. But this implies that S' \ {0} is totally real, and therefore the 
characteristic foliation has no singularities outside the origin. Thus, Theorem [1] holds under the 
assumption that is a generic real analytic diffeomorphism with D(p(0) symplectic. 

In the remaining part of this section we outline Bruno's algorithm, while the actual numerical 
calculations for (fTi|) and (p8|) are presented in the next section. 

Let X he a real analytic vector field on M? given by 

/ i = Ej+k>i(^jkt^s'' = tfi{t,s) 

We write 

f,{t,s) = Y,f3Qit^^)'^^ (30) 
Q 

where Q = {qi,q2) is the multi-index with integer entries, and {t,s)'^ = t^^s'^'^. The support D of 
X is the set of points Q = {qi, q2) in with integer coordinates such that |/iq| -|- |/2q| 7^ 0. The 
Newton polygon F is defined as the convex hull of the set 

P= U {(91,92) : 91 > 91, 92 > 92} • 

(gi,g2)eD 

The Newton diagram, or the open Newton polygon in the terminology of L of is the 
union of the compact edges of the Newton polygon T. This can also be obtained as follows. Let 

g2* = min{g2 : (91,92) G D} and qu = minjgi : (gi,g2*) G D}. The point := {qu,q2*) is the 
left boundary point of the intersection of D with the horizontal support line §2 = 92*- Consider 
the non- vertical support line L for D through F^'^'' which contains at least one other point of D. 
Denote by the left boundary point of the intersection of D with L. Consider now the support 
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line through r2'^^ which contains a point of D different from T^'^ . Continuing this procedure to 
the point Q* = {q^, g^g) which is the lowest point of D on the left vertical support line of D, i.e., 

= min{gi : ((71,(72) G D} and (^2 = min{g2 : {Qi,Q2) S D}. Denote this last point by r[.°\ 

For every j = 1, fc — 1 we denote by T^^'^ the edge joining the vertices T^J'^ and T^*^^- Thus by 

construction, the points vf'^ and T^^^ are joined by the Newton diagram f . 

(k) 

For each element F^- (j is enumeration as described above. A; = for vertices, k = 1 for 
edges) of the Newton diagram, there is a corresponding sector Uj in the phase space M^^ , so 
that together they form a neighbourhood of the origin (here boundaries of the sectors are not 
necessarily integral curves) . In each Uj one brings the system to a normal form, and in Uj one uses 
power transformations (quasihomogeneous blow-ups) to reduce the problem to study elementary 
singularities of the transformed system. After that the results in each sector can be glued together 
to obtain the overall phase portrait of the system near the origin. 

Several theorems on normal forms can be used to determine the phase portrait in a neighbour- 
hood of an elementary singular point. We state those of them that will be relevant to us. 

Principal Normal Form [Ul Ch. II, §1, Thm 2, p. 105] .• Given the system 

Xi = XiXi + o-jXj.i + ^ fiqX'^, i = 1, 2, 

|Q|>1 

with |Ai| + IA2I 7^ 0, there exists a change of coordinates Xi = yi + X]|q|>i ^iQ^^ > ^ = 1)2, such 
that in the new coordinates the system becomes 

where Ni = {qi > -l,qj > (j / i),qi + q2 > 0}, and giQ = unless {Q, A) = qiXi + q2X2 = 0. 

Consider the system of two differential equations in two variables of the form 

Xi = XiXi + Xi'^ figX'^ = XiXi + Xifi, i = l,2, A = (Ai,A2)/0, (31) 

QGV 

where V C is to be given. 

Second Normal Form [6l Ch. II, §2, Thm 1, p. 128]; Let R* and i?* be two vectors in 
contained in the second and the forth quadrant respectively, and let the convex cone V bounded by 
R* and R^ containing the first quadrant be such that its angle is less than n. Then system ([3T]) 
can be transformed by a formal change of variables 

Xi = yi + yi^h^QY^, i = l,2, (32) 

into a normal form 

yi = KVi + Vt'^giQY'^, i = l,2, (33) 

where Q £ Y, and giq = if {Q, A) ^ 0. 

Third Normal Form [6l Ch. II, §2, Thm 2, p. 134]; Suppose that the right-hand side of ([3T]) 
contains the series in integral, non-negative powers of X2, with the support entirely contained 
within the sector 

Y = {Q:Q = aiR* + 02^?*, ai, 02 > 0} , 
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where the vectors R* = {rl,r2) and = (ri^,,— 1) are such that rl < < r^, ri^ > 0, and 
ki/^2l < '"1*- This means that fig in fi{X) vanishes unless the vector Q lies in the sector 



Y = {Q:Q = aiR* + 02 ■ (1, 0), ai, 02 > 0} 



Denote by iV{X) the class of such series fi. The coefficient /2Q in f2{X) of (I3ip will vanish 
unless the vector Q lies either in iV, or along the ray {q2 = —1, qi > Denote the class of 

such series by 2V(X). Then the following holds: if the series fi in ()3ip are of class iV{X), then 
there exists a formal change of coordinates ()32p . where the hi are series of class iV(Y), which 
transforms (fSTj) into a system ([33l) in which the Qi are series of class iViY) consisting only of 
resonant terms gigY^ with {Q,A) = 0. 

In the cases which we will consider below, normalizing changes of coordinates will be analytic 
or at least C°°-smooth local diffeomorphisms, see [6j. This is sufficient in order to study a local 
topological behaviour of orbits. 

We now consider some important special cases. For the vertex Q = the corresponding 
sector is defined as 

(e) = {(t,s)GM2^|(t,s)|^*<e, < e, |i <e,\s\<e}, (34) 

where R* and i?* are the unit (i.e., their coordinates are coprime integer numbers) directional 
vectors of the edges adjacent to Q. The system can be brought to the Principal or the Second 
Normal Form by introducing a new time variable ri so that dri = {t, s)^dT. A particularly simple 

case occurs when Q = (gi, 92) = F^*^^ is the first or the last point of f and is not contained in the 
first quadrant. In this situation one of the coordinates of Q equals —1. Say, if ^2 = —1) then one 
takes iZ* = (0, 1), and the corresponding normal form has vertical integral curves. It follows that 
the original system ()29p in the sector 



U^{e) = {(t,s) G : |(t,s)|^* < e,\s\ < e} 

does not have any integral curves terminating at the origin. Similarly, if gi = —1, then R* = (1,0), 
and again in 

U*{e) = {{t, s) G : |(t, s)|^* < e, \t\ < e} 

the system does not have any characteristic orbits. 

Suppose now that F^-^'' is an edge of f . Let R = (ri,r2) be the unit directional vector of F^-^\ 
The corresponding sector in the phase space is given by 

Uj{e) = {{t,s) eM."^ : e <\{t,s)\^ <l/e, |t| < e, |s| < e}. (35) 
Consider the power transformation is given by yi = t^^ s^^ , 2/2 = t"^^ s^^ , where 

A=(^' (36) 
\ ri r2 J 

is a matrix with integer coefficients and determinant equal to 1. In the matrix form, we can write 

X = {t,s), 



Q 

Fn 



qi 
q2 

flq 
f2q 
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Then ()29p can be given by 

(InX) = FqXQ, (37) 

where X'^ = t^^s'^^. The power transformation can be expressed now as Y = X^ taking ([37]) into 

(iny) = F'q,YQ\ 
Q'eD' 

with Y = (1/1,2/2), Q' = {A^)^^Q, D' = (A*)^-'^D (the superscript t stands for transposition), 
and Fq, = AFq. After division by the maximal power of yi one obtains a new system. Here 
the 2/2-axis corresponds to {t = s = 0} in the original coordinates, and therefore one needs to 
investigate the new system in a neighbourhood of the 7/2-axis. The singularities of the new system 
are simpler than those of the original system, and therefore, an induction procedure can be used. 
Quite often the topological behaviour of the system in l^ji^) can be determined by considering 
the truncation of the system which is defined by taking the sum in (jSOp only over the vertices 
contained in T^}\ 



6. Phase portrait of the standard umbrella 

Since the standard umbrella corresponds to the non-generic case where (j) is the identity map, 
we study its characteristic foliation separately. We rewrite system ()14p in the form 



i = t(-3t2 - - 3t^) = tfi{t, s 
s = s(s2 + 4*2 + 7^4) = sf2it,s). 



and set 

Q 

where Q = (91, 92) is the multi-index with integer entries, and {t,s)^ = t'^^s'^'^. 

The Newton diagram f consists of two vertices r[^^ = (2, 0) and r2^^ = (0, 2) and the line 
segment (edge) T^^^ between them (see Fig. [T|). For each element of the Newton diagram (the two 
vertices and the edge), there is a corresponding sector in the phase space so that together 

they form a neighbourhood of the origin. Accordingly we consider 3 cases. 

Case 1. First consider the vertex (2,0). We define i?* = (1,0) and R* = (-1,1). The first 
vector arises since F^*^^ = (2, 0) belongs to the horizontal support line of D, and the second one is 
the unit vector in the direction of the edge in f starting from F^'^'' . We can make the change of 
time dri = t^dr. This yields the system 



The Newton diagram D corresponding to (j39p has vertices (—2,2) and (2,0), in particular, it is 
contained in the sector V (with the angle < vr) bounded by the rays generated by and R* . 
Therefore, for sufficiently small e, in the sector 

Uf^ = {(t, s) G m2 : (|t|, |s|)«* < e, (|t|, |s|)^* < e} = {\t\ < e, \s\ < e\t\}, 
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Figure 1. The Newton diagram for ([SS]) . 



there exists a smooth change of variables {t, s) putting the initial system to the Second Normal 
Form of Bruno. In the new coordinates the system has the form 

m = -3yi + yi E giQiviiVi)^ .^qn 
y2 = 4y2 + y2 E52Q(yi,y2)'^, 

where the coefficients giq and g2Q are all zero except those for which — 3gi + 4(72 = 0. The line 
L := {— 3yi + 4?/2 = 0} determined by the linear part of system (|40p intersects the interior of the 
sector V (see Fig. [2|). It follows (see Bruno p. 132) that the system defined by (|10]) . and hence 
by (f39|) . is a saddle, i.e., each ray {iy^ > 0,j = 1,2} is an integral curve, and in each quadrant 
in R^, the integral lines are homeomorphic to hyperbolas. This is the description of system (|17p 
in sector uf'* . 

Case 2. Consider now the second vertex (0,2). Again, following [6] define i?* = (1,-1) and 
R* = (0, 1). The first vector is the direction of the edge starting at the vertex (0, 2); the second 
one corresponds to the vertical support line and is imposed since this vertex is the left boundary 
point of the Newton diagram. The corresponding sector where the change of dependent variables 
will be performed is given by 



U^^^ = {{t,s) : (|t|,|s|)^* <e, {\t\,\s\)^* < e} = \ \s\ < e, \s\ > 
The change of time dri = s'^dr transforms system (I17p into 

] |i = s + s(4t2s-2 + 7t4s-2). 



} 



(41) 
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Figure 2. Case 1 for ([Ml)- 



As above, there exists a smooth change of variables (t, s) putting this system to the second normal 
form: 

yi = -yi + yiT,9iQ{yi,y2)'^ 
m = y2 + y2E52Q(yi,y2)'^, 

where the coefficients qiq and g2Q are all zero except those which belong to the line L := {—qi + 
(72 = 0}. This line intersects the sector V bounded by i?* and R* which implies that this system 
is again a saddle. This gives the phase portrait of (fT7|) in sector Wg^^ 



u[^^ , which is the complement of L(['^^ U ^2"'' • We make the following change of variables 



Case 3. The remaining case of the edge between (2, 0) and (0, 2) will correspond to the sector 

. We m 

= lu (42) 

y2 = t s- 

In the matrix form, we write X = {t,s), and the change of variables (j42p can be expressed as 
Y = with the matrix of exponents 



A 

Then system (jl7p takes the form 



1 
-1 1 



yi = yi{-3yf - yfyi - 3yf) 
y2 = y2{7yl + 2ylyj + Wyf) 




Figure 3. Case 2 for §8^. 



The edge of f becomes vertical in the new system. Performing as above a change of time, we may 
divide both sides by yf to obtain 

/ m = -3yi - yiyi - 3yf = yi{-3 - yl - 3yj) , . 

1 m = 7y2 + 2yl + IQyly^ = y2(7 + 2yl + lOy?). ^ 

Under the change of variables (j42p . the line yi = corresponds to the origin, and therefore, we 
are interested in the integral curves of system (j43p that intersect the line yi = at points with 
2/2 7^ 0. The set {yi = 0, ±y2 > 0} are integral curves of ([MD, but they correspond to t = s = 
in the original system. According to Bruno ([6], p. 141), the points on the y2 axis can be either 
simple points, in which case the integral curves of (j43|) near such points are parallel to the y2-axis, 
or singular points. The truncation of system (I43p (see the end of the previous section) contains 
only the terms that correspond to the edge under consideration and its vertices, and thus has the 
form 

yi = yif},{y2) ^^^^ 
m = ^2/20 (2/2), 

where f'^^{y2) = 7 + 2y| (we follow the notation of [6j). Singular points are determined from the 

equation f'^^{y2) = 0. In our case f'^^{y2) is strictly positive. Therefore, in all points with 

yi = 0, y2 7^ are simple points. From this we conclude that in the sector no integral curves 
of system ()17p intersect the origin. 

With this information the integral curves in all sectors can be glued together. It is readily 
verified that the phase portrait of system (I17p is in fact a saddle, the integral curves in each 




Figure 4. Phase portrait of (f38|) . 



quadrant of are homeomorphic to hyperbolas and do not intersect the coordinate axes (see 
Fig.H. 

7. Phase portrait of umbrella in general position 

We now perform similar calculations for the algorithm to determine the topological structure 
near the origin of the dynamical system defined by ()28p . First of all we represent it in the canonical 
form 

J i = t{-2gi2S + ao2t"^s2 - 8522*^ + •.•) 

\ s = s{4gut^ + Puts + /3o3s' + 6512*^5"^ + •.•)• ^ ' 

The Newton diagram F consists of 3 vertices (—1, 2), (0, 1) and (4, —1), and the two edges between 
them (Fig. E]). Five cases should be considered each corresponding to a vertex or an edge of F. 

Case 1. Vertex (4, —1). This is described in ^ p. 138] as Type L As as result, we immediately 
obtain the behaviour of integral curves of the system. Namely, in the sector 

= {(t,s) eM2 . (|t|,|s|){i.o) (|^|j^|)(~2,i) <e] = {\t\ <e, \s\ < e\tf} 
the integral curves are vertical, in particular, there are no curves terminating at the origin. 
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Case 2. Vertex (—1,2). Again the same analysis works here. Since (—1,2) is the end point of 
r, i.e., of Type I in [(3, p. 138], it foUows that in 

Up = {{t, s) G m2 . (1^1, |s|)(0-i) < e, (|t|, |s|)(^'~2) <e] = {\s\ < e, \t\ < e\s\^} 

the integral curves are horizontal, and no curves terminate at the origin. 

Case 3. Vertex (0,1). This is Type III in [6l p. 139], p. 123. After a change of time so that 
dri = sdr, the system takes the form 

i = t{-2gi2 + ao2t~^s - 'ig-22t^s~^ + ...) . . 

s = s(45iit's~' + Pi2t + /3o3S + 6(7i2t^s-2 + ...). ^ 

There are two sectors which can be assigned to vertex (0,1). One of them is determined by 
i?* = (2, —1) and R* = (—1, 1), and equals 

= {{t, s) G ]R2 : (|t|, l^l)^* < e, (|t|, |s|)«* < e]. 

We may apply here the Second Normal Form of Bruno. Since we consider a generic case, we have 
Ai = — 2(7i2 7^ 0. Further, A2 = 0, because the second equation has no free term. Recall that we 
use the notation A = (Ai, A2). The line L determined by 

L = {Q = {qi,q2) e M2 :< Q, A >= 0} = {q^ = 0} (47) 
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enters the interior of the sector bounded by R* and i?*. It follows that in U2^^ there are no integral 
curves terminating at the origin. 

On the other hand, we may use the Third Normal Form of Bruno for (|46p . It is valid on a 
bigger domain, namely, on 

= {{t,s) gM^ : (|t|,|s|)-^* <e, |s| < e} = {\t\'^ <e|s|,|s| < e}- 

The region of the (t, s)-space where the dynamics takes place is given by 

2V = {Q : Q = aiiJ* + 02 • (0, 1), ai,a2>0}. 

Now the line L determined from ()47p enters 2V along its boundary, the s-axis. In general, this 
yields a complicated behaviour of the system in 2^2^^ ■ In fact, there are four possibilities as 
described in [Gj p. 134 Case c)]. So which case is it? The salvation comes from Case 2 above: 

it describes the behaviour of the system in U^^^ (which is a subset of 2^2^^^ and a neighbourhood 
of the s-axis). According to Case 2, the integral curves are horizontal near the s-axis, which 

eliminates all possibilities but one. We conclude that no integral curves enter the origin in 2^2^^ ■ 
Case 4- Edge connecting (0, 1) and (—1,2). The corresponding sector is defined by 

Z^W = {QeR'^:e< (|t|, |s|)(-^'i) < e'^} 

(see [Sl p. 139]). This case is subsumed by Case 3 above because Z^2^^ C 2^2^^^ in a suitable 
neighbourhood of the origin. 

Case 5. Edge connecting (0,1) and (4,-1). We will consider the truncation of system (f28]l . 
i.e., we keep only terms that are related to the edge under consideration. We have 

i = t{-2gi2S - 3522*^) z.^x 
s = s{-Agnt^ + 6gi2t^s-^). ^^^^ 

The directional vector is R = (—2, 1), and the sector in which the dynamics should be understood 



IS 



U[^'^ = {{t,s) ■.e< (t,s)(~2'^) < <e} = {e|tp < \s\, \s\ < \t\ < e, \s\ <e}. (49) 



We need to make the following change of coordinates: 

yi = t 
y2 = t-^s, 



(50) 



which corresponds to the matrix 



A 



1 
-2 1 



In the new coordinates system (|48p becomes 

yi = yi(-25i2y?y2 - 35222/?) 
m = y2{^9i2yly2 + (6522 + 4511)2/1 + 69122/12/^^). 

We divide by the maximal power of yi, which equals 2 in this case, by performing the change of 
the independent variable: dri = yldT. This yields 

m = 2/1 (-2^12^2 - 3522) /g2) 
2/2 = 2/2(45122/2 + (6522 + 4c/ii) + 65122/2"^). 
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Figure 6. Phase portrait in y-coordinates for gi2 > 0. 



This is the system of Type I in [6], p. 125]. The y2-axis is an integral curve, but it corresponds to 
the origin in PH]) . Near points where i{^,y2) = 4(712^2 + (6522 + '^9ii)y2 + 6(712 is not zero, the 
integral curves are parahel to the y2-axis, so these are the curves that do not terminate at zero 
in the original system. We need to investigate the situation near points where ^(0)2/2) = 0. For 
this we solve the quadratic equation 

2^12^1 + (3522 + 2511)2/2 + 3c/i2 = 0. (53) 

The discriminant of this equation is 

V = Agl^ + 9gi2 + 12511522 - 245^2 • 

Since Agl^ + 9^22 > 125111522, it follows that V > 24511522 - 245^2 = 24A > 0. Here A is defined 
by dSI). Thus, equation (j53|) always has two simple roots: 

-(3g22 + 2gii) ± ^/Agf^ + 9g^2 + 12giig22 - 24ff^2 

4c/i2 



(since we consider the generic case, we can assume that gi2 7^ 0). We point out that c± are either 
both positive or both negative. 
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We need to investigate the dynamics near each point c±. For that we first need to translate c± 
to the origin via 

zi = yi, y2 = c± + Z2. 
In the new coordinates the system becomes 

f ii = zi{-{2gi2C± + 3522) - '2gi2Z2) /g^N 

\ Z2 = Z2{{8gi2C± + 6522 + 4511) + 4:gi2Z2)- 

This is a system for which the origin in an elementary singularity (the linear part is not zero) . To 
determine the dynamics we need to understand the sign of the coefficients of the linear part, i.e., 
of 



and 




Figure 7. Phase portrait of gu > 0. 
Claim. Ai and A2 are of the opposite sign both for c+ and c_ . 

First note that Ai and A2 depend only on the coefficients gjk, i.e., only on the linear part of 
the map tp o (f). Therefore, it is enough to prove the claim for linear symplectomorphisms. If (j) is 
the identity map, then it is easy to see that Ai and A2 are of the opposite sign. 
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Suppose that for some hnear symplectic map (po, the sign of Ai and A2 is the same. Since the 
symplectic group is connected, there is a path 7 C Sp(4, M) connecting the identity and (po, and 
since Xj depend continuously on <f>, there exists a symplectic map on 7 for which one of the A^ is 
zero. Since P > 0, it has to be Ai. So —^922 + dn = ±^\/P. Therefore, 

4gj^ - 12511922 + Qgh = ^ali + 9^22 + 12(7ii522 - 245^2. 
This implies that A = - contradiction. This proves the claim. 

Since Xj are of different sign, it follows that both for 0+ and c_, system (|54p is a saddle at 

the origin. Now we are able to describe the overall dynamics in In (yi,?/2)- coordinates we 

have the following: y2-axis as well as the lines y2 = C4. and y2 = c_ are the integral curves. More 
precisely, the integral curves are six half-lines: Li = {(yi,c+),yi > 0}, L2 = {(yi,c+),yi < 0} , 
L3 = {(yi,c_),yi > 0}, L4 = {(yi,c_),yi < 0}, L5 = {(0,7/2) : y2 > c+}, Lq = {(0,7/2) : y2 < C-}, 
and one interval / = {(0,1/2) : min{c_,c+} < 1/2 < max{ c_,c+}} . The phase portraits near the 
points (0,c+) and (0, c_) are saddles, whose orbits in between the lines 2/2 = c+ and 1/2 = c_ are 
glued together, and are asymptotic to Li , L3 or to L2, L4 ; they do not touch / . Other orbits 
are asymptotic to L2, or to L5, Li or to Lg, -L4 or, finally, to Lg, -L3 (see Fig. [6|). Going back 
to the original system via the inverse transformation to ()50p . we see that the y2-axis blows down 
to a point, and we have four integral curves s = c±t^ entering the origin, while other integral 
curves are contained in the compliment of these two curves. Now, if we choose e > sufficiently 
small in ()49p . we see that both curves s = c±t^ enter IA^\ This completes Case 5. 

Now if we combine all 5 cases together, and glue the integral curves from all cases, we see that 
the phase portrait at the origin of system (I28p is a saddle (Fig. [7]). With this analysis we can now 
conclude the proof of Proposition [T] Indeed, let 71 and 72 be the curves s = c±t^. If if is a small 
compact not contained in the union of 71 and 72, then one of the hyperbolas of the characteristic 
foliation will touch K at some point. This proves Proposition [T] 
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